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Abstract 

We review the Floquet theory of Unear differential equations with periodic coef- 
ficients and discuss its applications to neutrino oscillations in matter of periodically 
varying density. In particular, we consider parametric resonance in neutrino oscil- 
lations which can occur in such media, and discuss implications for oscillations of 
neutrinos traversing the earth and passing through the earth's core. 
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1 Introduction 



All oscillating systems are very much alike, and there are many similarities between oscil- 
lating neutrinos and, e.g., pendulums or electromagnetic circuits. In particular, neutrino 
oscillations in vacuum or in matter of constant density are analogous to oscillations of a 
simple pendulum; resonantly enhanced neutrino oscillations in a matter of monotonically 
varying density (the Mikheyev-Smirnov-Wolfenstein (MSW) effect [|l], 0) are similar to os- 
cillations of two weakly coupled pendulums of slowly changing lengths P, §]. It is therefore 
natural to ask the following question: Are there any other resonance phenomena in mechan- 
ics or electromagnetism that might have analogues in neutrino physics? 

One such phenomenon is parametric resonance. The parametric resonance can occur in 
dynamical systems with time- varying parameters when there is a certain correlation between 
these variations and the values of the parameters themselves. Best studied is the parametric 
resonance in systems with periodically varying parameters. While the periodicity makes it 
easier to satisfy the resonance conditions and also simplifies the analysis, it is not really 
necessary: parametric resonance can occur even in stochastic systems (see, e.g., 0). In the 
present paper we will concentrate on the systems with periodically varying parameters. 

A textbook example of a system in which the parametric resonance can occur is a pen- 
dulum with vertically oscillating point of support 0. Under certain conditions topmost, 
normally unstable, equilibrium point becomes stable. The pendulum can oscillate around 
this point in the upside-down position. Another example, familiar to everybody, is a swing, 
which is just a pendulum with periodically changing effective length. It is the parametric 
resonance that makes it possible to rock a swing. 

What would be an analogue of the parametric resonance for neutrino systems? Since 
matter affects neutrino oscillations, periodically varying conditions can be achieved if a 
beam of oscillating neutrinos propagates through a medium with periodically modulated 
density. For certain relations between the period and amplitude of density modulation 
and neutrino oscillation length and mixing angle, the parametric resonance occurs, and 
the oscillations can be strongly enhanced. The probability of neutrino transition from one 
flavor state to another may become equal to unity. This phenomenon is very different from 
the MSW effect. Indeed, at the MSW resonance the neutrino mixing in matter becomes 
maximal {6m = 7i"/4) even if the vacuum mixing angle 6o is small. This leads to large- 
amplitude neutrino oscillations in a matter of constant density equal (or almost equal) to 
the resonance density, or to a strong flavor conversion in the case of matter density slowly 
varying along the neutrino path and passing through the resonance value. 

The situation is quite different in the case of the parametric resonance. The mixing angle 
in matter does not in general become large (there is no level crossing). What happens is 
an amplification of the transition probability because of specific phase relationships. Thus, 
in the case of the parametric resonance it is the phase of oscillations (rather than their 
amplitude) that undergoes important modification. The total flavour conversion can take 
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place even if the mixing angles both in vacuum and in matter are small. 

The possibility of the parametric resonance of neutrino oscillations was suggested inde- 
pendently in and . In these papers approximate solutions for sinusoidal matter density 
profile were found. In also an exact analytic solution for the periodic step-function 
( "castle wall" ) density profile was obtained. Parametric effects in neutrino oscillations were 
further studied in where combined action of the parametric and MSW resonances and 
possible consequences for solar and supernova neutrinos were considered. In this paper also 
the stochastic parametric resonance in neutrino oscillations was briefiy discussed. 

Although the parametric resonance in neutrino oscillations is certainly an interesting 
physical phenomenon, it requires that very special conditions be satisfied. Unfortunately, 
these conditions cannot be created in the laboratory because this would require either too 
long a baseline or neutrino propagation in a matter of too high a density. Until recently it 
was also unclear whether a natural object exists where these conditions can be satisfied for 
any known source of neutrinos. This situation has changed with a very important observa- 
tion by Liu and Smirnov |Tl| (see also who have shown that the parametric resonance 



conditions can be approximately satisfied for the oscillations of atmospheric z/^ into sterile 
neutrinos Us inside the earth. The density profile along the trajectories of neutrinos crossing 
the earth and passing through its core (mantle-core-mantle) is to a good approximation a 
piece of the periodic step-function profile, and their oscillations can be parametrically en- 
hanced. Even though the neutrinos pass only through "1.5 periods" of density modulations 
(this would be exactly one period and a half if the distances neutrinos travel in the mantle 
and in the core were equal), the parametric effects on neutrino oscillations in the earth can 
be quite strong. Subsequently, it has been pointed out in |jl3| that the parametric resonance 
conditions can also be satisfied (and to even a better accuracy) for the 1^2 z/g oscillations 
in the earth in the case of the z/g - i^^(t) mixing. This, in particular, may have important 
implications for the solar neutrino problem. The parametric resonance in the oscillations 
of solar and atmospheric neutrinos in the earth was further explored in a number of papers 

g ig ig 0, 0, [T|,g. 



In the present paper we review the Floquet theory of linear differential equations with 
periodic coefficients and consider its applications to neutrino oscillations and, in particular, 
to oscillations of neutrinos inside the earth. The paper is organized as follows. In sec. 2 
we briefly review the Floquet theory and its applications to the analyses of the stability 
of the solutions. In sec. 3 we discuss the peculiarities of the Floquet theory in the case 
of the time dependent Schrodinger equations with periodic Hamiltonians. In sec. 4 we 
consider applications of the Floquet theory to neutrino oscillations in matter of periodic 
step function ("castle wall") density proflle. In sec. 5 we review the implications of the 
parametric resonance of neutrino oscillations for neutrinos traversing the earth and passing 
through its core. In the last section the results are discussed and the conclusions are given. 
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2 Differential equations with periodic coefficients 



We shall now briefly review the Floquet theory of systems of linear differential equations 
with periodic coefficients. More detailed discussion can be found, e.g., in [0, pT[ ]. 

2.1 Preliminaries 

Let us start with a few well known fact from the general theory of linear differential equa- 
tions. Consider a system of n homogeneous linear differential equations 

^ = A{t)i,, (1) 

where is an ?7,-component column vector, ip = {ipi, ...jipn)'^ , A{t) is an tt, x n matrix 
with piecewise continuous elements, and overdot denotes differentiation with respect to t. 
Eq. dll) has n linearly independent continuous nontrivial solutions ip^^^t), j = 1, n. From 
linearity of (|l]) it follows that any linear combination of the solutions is also a solution. Any 
set of n linearly independent solutions il)'^^\t) of (|I]) forms the so-called fundamental set, 
and a matrix whose columns are i})^^^ is called the fundamental matrix. Given an initial 
condition ip{tQ) = ipQ, eq. (0) has a unique solution ip(t). Any solution of eq. (|^) can 
be represented as a linear combination (with constant coefficients) of the solutions forming 
a fundamental set, or equivalently as a product of a fundamental matrix and a constant 
vector. 

Let ip(t) be the solution of ([|) with the initial condition ip(to) = ipo- Let us introduce 
the evolution matrix U{t, to) through the relation 

m = Uit,to)^Po. (2) 

From the definition of the evolution matrix it immediately follows that 

U{t, to) = U{t, ti)f/(ti, to) , U{to, to) = / , (3) 

where I is the nxn unit matrix. It is easy to check that the columns of U (t, to) are solutions 
of eq. (|l|) with the initial conditions ^p)(0) = 5ij. Thus f/(t,to) is a fundamental matrix, 
and any solution of (|1|) can be written in the form (Q). The determinant of U{t,to) is given 
by 

det[U{t, to)] = exp |^ tr[A{t')] rft'j . (4) 

Eq. @ follows from two facts: (1) the derivative of a determinant is the sum of n determi- 
nants formed by replacing elements of one row of the original matrix by their derivatives, and 
(2) the columns of U (t, to) are solutions of eq. (|l]). Since the elements of A are non-singular, 
det[U] does not vanish, i.e. the matrix U{t,to) is non-singular. From (|^) one finds 

U{h,h)-' = U{h,t2). (5) 
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Without loss of generality, one can always choose the matrix A{t) to be traceless. Indeed, 
substituting 

^(t) = e"(*V W , a{t) = ^f tr[A{t')] dt' (6) 
one finds that tp' satisfies 

^' = A'it)^' , A'it) = Ait) - -trlAit)] , (7) 

n 

i.e. A' is traceless. In what follows we will be always assuming that the transformation (P) 
has been performed, i.e. will be considering only traceless matrices A. Eq. (^) then yields 

det[U{t,to)] = l. (8) 

Notice that in general U{t, to) is not unitary. 

2.2 Floquet theory 

Let us now turn to the case of differential equations ([l|) with periodic coefficients, 

A{t + T)= A{t) , (9) 

where T is the period. Hereafter we shall be always assuming that (Q) is satisfied, without 
specifying this each time explicitly. From the periodicity of A{t) it follows that if iplt) 
is a solution of (|1|), so is iplt + T). Consider the solution iplt) with the initial condition 
ip{s) = ipQ. We have tpit + T) = U{t + T, s)ipo- On the other hand, since ip(t + T) is also a 
solution at t, we have ip(t + T) = U {t, to)i/j(to + T) = U(t, to)U(to + T, s)iIjo. Equating these 
two expressions for ijj{t + T) one obtains 

U{t + T,s) = U{t,to)U{to + T,s). (10) 

This is a very important property of the evolution matrix of differential equations with 
periodic coefficients. In particular, taking s = to + T and s = to = one obtains from ([TOD , 
respectively, 

[/(t + T,to + T) = f/(t,to), (11) 
f/(t + T, 0) = f/(t, 0) f/(T, 0) . (12) 

The first of these equations means that the evolution matrix does not change if both its 
arguments are shifted by the period T (and, by induction, by any integer number k of 
periods). From the second equality it follows, in particular, that the matrix of evolution 
over k periods satisfies 

f/(fcT, 0) = f/(T, 0)^ (13) 

The matrix of evolution over one period U{T, 0) = Ut plays a very important role in the 
theory of differential equations with periodic coefficients; it is called the monodromy matrix. 
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As has already been pointed out, if ip(t) is a solution of ([l|), so is ipit + T). This does 
not in general mean that ilj{t + T) = i.e. the solutions of the equations with periodic 
coefficients are not in general periodic. There are, however, solutions which satisfy 

^{t + T) = am, (14) 

i.e. they are multiplied by a number when t is shifted by the period. Such solutions are 
called normal; they play an important role in the analysis of the stability of the solutions 
ofeq. (I). 

To analyse the properties of the normal solutions, let us show that the evolution matrix 
for a system of linear differential equations (0) with periodic coefficients can be written as a 
product of a periodic matrix and an exponential matrix. Since the monodromy matrix Ut 
is non-singular, it can be represented as an exponential of another matrix: 

Ut = U{T, 0) = e^^ . (15) 

Let us now show that the matrix P{t) = ?7(t,0)e~^* is periodic. We have P{t + T) = 
U{t + T, 0)e-^(*+^) = U{t, 0)t/(T, 0)e-^^e-^* = U{t, 0)6"^* = P(t), where we have used (0) 
and ([TSl) . Thus the evolution matrix can be written as 

f/(t, 0) = P(t)e^* , P{t + T) = P{t). (16) 

From this expression it follows, in particular, that the vector x introduced through ijj = 
P{t)X satisfies the differential equation with constant coefficients: 

X = Bx- (17) 



Let 0Q be a (constant) eigenvector of the matrix B with the eigenvalue ajT ^. It is 
then also an eigenvector of Ut with the eigenvalue cTj = e"^ : 

= f rf^ ; ^T0i^'^ = = e-rf) . (18) 

The numbers (Tj are called the characteristic numbers and aj, the characteristic exponents. 
From (|1^), (|16D and ( [T8| ) it follows that any normal solution (f)^j\t) of eq. (|I]) can be written 

as 

^(^■)(t) = f/(t, O)0i^') = P(t)e"^(*/^VS'^ = P{t)(Tf^^i'^ ■ (19) 

Thus eigenvectors of Ut give rise to normal solutions. The monodromy matrix Ut has n 
eigenvalues. As follows from (H), they satisfy 

n 

n^. = i- (20) 

If all the eigenvalues aj are different, Ut has n linearly independent eigenvectors and there- 
fore there are n linearly independent normal solutions (j)^^\t) of system (|T]). Hence, they 
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form a fundamental set, and any solution of (|T|) can be written as a linear combination of 
the normal solutions (f)'^^\t) with constant coefficients. If Ut has repeated eigenvalues, the 
situation is more complicated and will be discussed below. 

It follows from eq. (0) that the characteristic exponents (or characteristic numbers) 
determine the boundedness of the normal solutions, and therefore of the general solution 
of eq. (|l|). Consider, for example, the case n = 2. Assume first that the characteristic 
exponents are real. The characteristic numbers ai and (T2 are then real, too. Eq. ( [2D| ) 
gives aia2 = 1, and if ai 7^ (T2 (i-e. they are simple eigenvalues of Ut), the absolute value 
of one of them is greater than one. It then follows from eq. (|T9|) that the corresponding 
normal solution is unbounded. If the characteristic exponents are purely imaginary, the 
characteristic numbers are of modulus one, i.e. ai = cr2, |c"i| = 1- In this case both normal 
solutions are bounded (the matrix P(t) in (|19|), being continuous and periodic, is obviously 
bounded). If the characteristic exponents are complex, the normal solutions are bounded 
when the real parts of all the characteristic exponents are non-positive, while if at least 
one of the characteristic exponents has a positive real part, there are unbounded normal 
solutions. This property holds in general, i.e. for an arbitrary n. 

Let us now discuss the case of repeated eigenvalues of Ut- Consider again the case n = 2 
as an example. Assume that the eigenvalues of Ut coincide, ai = a2 = u- The matrix Ut 
in this case has the general form Q 

It can be readily seen that for a = the monodromy matrix Ut has two linearly independent 
and orthogonal eigenvectors corresponding to the same eigenvlue a; they can be taken, e.g., 
as 00^^ = (0, 1)-^ and ^q^"* = (1,0)-^. However, for a 7^ there is only one eigenvector, ^q^"*. 
It gives rise to the normal solution (j)'^^\t) through eq. ([Tol) . This normal solution can be 
used as one of the basis solutions ip^^\t) constituting a fundamental set, ip^^\t) = (p^^^t). 
Let il)^'^\t) be another solution of (P, linearly independent from il}'^^\t). Then ip'^^^t) and 
ip'^'^\t) form a fundamental set. Since ijj^^\t + T) are also solutions, they can be written as 
linear combinations of ■ijj^^\t) and il)'^'^\t): 

^(2) (t + T) = aV^i) (t) + 6^(2) ^ ^22) 

Using the freedom of normalization of the vectors (f)Q \ one can always choose a' = a. Since 
0" is a double root of the characteristic equation of the monodromy matrix Ut, b' = cr. 

^Another possibility would be to have the zero element {Ut)21 instead of {Ut)i2, but the corresponding 
matrix can be reduced to that in (|2^) by a renumbering of the basis states. 

^An easy way to see this is to consider eqs. ( ^2| ) a.t t — and use the explicit form of the matrix 
Ut- 
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Eqs. (^) can therefore be rewritten as 

V^(i)(t + r) =a^«(t), (23) 

^(2) (t + T) = aij^^^ (t) + ai,^^^ {t) . (24) 

Let us introduce a matrix M which relates tpl^^ and ■i/'o^'' = ^q^"*: ipl^^ = Mip^^ ^ One can 
now find the relation between ip^'^\t) and ip^^\t): 

■ilj(^){t) = W{t)ilj'-^\t) , W{t) = U{t, 0)MU{t, 0)-^ . (25) 

Therefore ifj'^^^ (t + T) =W{t + T)^)^^'^ (t + T) = W{t + T)ai)^^^ (t) . On the other hand, from 
dH) and (HD, ^^^\t + T) = a^/'(i)(t) + a^p^^^t) = [a + (Tiy(t)]^/'(i)(t). Equating these two 
expressions for il)^'^\t + T) one finds 

(t + T) = W^(t) + - . (26) 
a 

Therefore the matrix F{t) defined through 

F{t) = W{t) --^ (27) 

is periodic with the period T. From eqs. ( [25| ) and ( pTj) we then find 

"a t 



cr T 



^^^^(t), F(t + T) = F(t). (28) 



This is the result we were looking for. Together with eq. ( ^31) it states that in the case of the 
double roots of the characteristic equation of Ut (i.e. in the case of coinciding characteristic 
exponents), a fundamental set can be chosen to consist of a normal solution and a solution 
which is a linear combination of a linearly growing and periodic functions multiplied by the 
normal solution. Thus, if a 7^ 0, there are unbounded solutions. As follows from (^), if 
a = 0, then tp^'^\t), being a product of a periodic matrix and a normal solution, is also a 
normal solution. Thus in this case both solutions forming a fundamental set can be chosen 
to be normal. This is in accord with the fact that for a = the matrix Ut in eq. (l2l|) has 
two linearly independent eigenvectors. 

It is instructive to see how the linear growth of ip^'^\t) arises from the general expression 
7/>(2)(t) = Pii^e'^^ipfK From (0) and (|l]) one finds 

A2 

-BT = Inf/T = ln[a(J + A)] =lna + ln(/ + A) = lna + A + ... , (29) 

2 3 

where 

A = 




^The matrix M is not uniquely defined, but this is unimportant for our purposes. 
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Since = 0, eq. (|9D gives ST = In a + A, hence 

Notice that the matrix A annihilates ^q^'*, and therefore ( ^OD does not contradict ip'^^\t) 
being a normal solution. Using eq. ( pO| ) one can find a simple representation of the matrix 
F{t) entering into eq. (^), 

F{t) = P{t)MP{t)-\ (31) 

from which the periodicity of F{t) is obvious. 

The above result can be generalized to the case n > 2. If the characteristic equation 
of the monodromy matrix Ut has repeated roots (i.e. some of the characteristic numbers 
coincide), a fundamental set can be chosen to consist of normal solutions and solutions 
which are linear combinations of polynomials in t and periodic matrices multiplied by normal 
solutions. 

We have seen that in the case when some of the characteristic exponents have positive 
real parts, there are exponentially growing solutions, while if some of the characteristic 
exponents coincide, there are in general polynomially growing solutions. The existence of 
such unbounded solutions signifies instabilities due to the parametric resonance. 

It should be noticed that instead of considering a system of n first order linear equations 
with periodic coefficients (0) one could equivalently consider one equation of the nth order. 
In particular, in the case n = 2, a general second order equation with periodic coefficients is 
obtained, which is called the Hill equation. A pendulum with vertically oscillating point of 
support mentioned in the Introduction is described by this equation. If the oscillations of 
the point of support are harmonic, the pendulum is described by the well known Mathieu 
equation. In the limit of small- amplitude oscillations of the point of support, the instability 
condition (the condition of exponential growth of the deviation from the equilibrium) is 

n.f^f, (3.) 

where u is the frequency of the oscillations of the pendulum in the absence of the motion of its 
point of support, and k is an integer. Eq. ( p2|) relates the frequency of the oscillations of the 
point of support Q at which the parametric resonance occurs to the oscillator eigenfrequency. 
In general, when the amplitude of the oscillations of the point of support is not small, the 
parametric resonance condition depends not only on the frequency of these oscillations but 
also on their amplitude. In this case there are resonance regions of parameters rather than 
resonance values 0, 21 . 



In real physical systems all parameters are, of course, finite; unboundedness of certain 
solutions of eq. (P with periodic coefficients just refiects the fact that in general the dynam- 
ics of real systems is only approximately described by linear equations. For large deviations 
from equilibrium, nonlinear effects become important and eqs. (|1|) have to be modified. 
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There are, however, cases when the solutions are always bounded even in the linear regime, 
and the description by linear equations can be exact. Nevertheless, parametric resonance is 
possible in such systems as well. One example of such a situation is given by Schrodinger 
equations with periodic Hamiltonians, which we discuss next. 

3 Schrodinger equations with periodic Hamiltonians 

If the matrix A{t) in eq. (|l|) is anti-Hermitian, the system of equations (|lD, (§) can be 
written as a Schrodinger equation with a periodic Hermitian Hamiltonian: 

i^ = H{t)tij, n{ty = n{t), n{t + T) = n{t), (33) 

where Hit) = iA(t). In the case of constant H, eq. (^) describes oscillations between the 
components of ip characterized by n— 1, in general different, frequencies (out of n eigenvalues 
of Ti only n — 1 are independent since trTi = 0). In particular, spin precession in a constant 
magnetic field or neutrino oscillations in vacuum or in matter of constant density are de- 
scribed by such an equation. Eq. (|33D with time-dependent periodic Hamiltonians describes 
many physical systems, e.g., atoms in a laser field or electron paramagnetic resonance. It 
also describes neutrino oscillations in a medium of periodically modulated density. 

Because of the Hermiticity of TC(t), the evolution matrix U{t,to) of eq. ( ^31) is unitary, 
i.e. the norm of the vector is conserved. Therefore all the solutions of eq. ( P^ ) are 
bounded. The parametric resonance in the systems described by eq. (|33D has therefore some 
peculiarities, which we shall discuss below. 

One consequence of the unitarity of the evolution matrix U{t, to) is that all the char- 
acteristic exponents are purely imaginary. In addition, since the polynomial growth of the 
solutions is not allowed in this case, even in the case of repeated roots of the characteristic 
equations of the monodromy matrix Ut there are n linearly independent normal solutions 
which form a fundamental set. This actually directly follows from the fact that a unitary 
n xn matrix has exactly n linearly independent eigenvectors, irrespective of whether or not 
all the roots of its characteristic equations are simple. 

We shall now discuss general properties of the solutions of eq. (^), using again the case 
n = 2 as an example. The results, in particular, will apply to the problem of two-flavour 
neutrino oscillations in a matter of periodically modulated density. 

Let us start with a few general remarks about the solutions of the Shrodinger equation 
with time dependent (but not necessarily periodic) Hamiltonian in the case n = 2. First, 
we notice that without loss of generality the Hamiltonian Ti.{t) can be considered to be real. 
Indeed, in the case of a complex Hamiltonian, a rephasing of the components ?/'i_2 of by 
the factors exp(±i/5(t)/2) where /5(t) = arg[Hi2{t)] transforms the Hamiltonian to the form 

- ( I! ) (-) 
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with real A{t) and B{t). Notice that this rephasing transformation preserves the trace of 
the Hamiltonian. 

Next, we notice that the Hamiltonian 7i{t) in can be written as 

n{t) = B{t)ai - A{t)a^ , (35) 

where (Tj are the Pauli matrices. Thus 'H{t) anticommutes with a2- From this fact it 
immediately follows that if il){t) = (^/'i(t), ■?/'2(t))"^ is a solution of the Schrodinger equation, so 
is ipit) = —ia2ip{t)* = {—■ip2(t)*,ipi{ty)'^. It is easy to see that ipit) and ipit) are orthogonal 
and therefore linearly independent. Thus, if one nontrivial solution of the Schrodinger 
equation in the case n = 2 is known, it automatically gives another nontrivial solution, 
linearly independent from the original one. The solutions ip and ip form a fundamental set 
and therefore knowledge of a single nontrivial solution of the Schrodinger equation allows 
one to obtain the general solution. 

We now turn to the study of the Schrodinger equation with periodic coefficients in the 
case n = 2. Since the monodromy matrix Ut = exp{BT) is a unitary 2x2 matrix, it can 
be written as 

Ut = Y - iaX = exp[-i(crX)$] (36) 
with real parameters X and Y (or X and $) Q which satisfy 

Y^ + X^ = l, cos<l> = F, sin$ = |X|; X = ^. (37) 

I I 

It should be noted that the form of the monodromy matrix in eq. ( pBf ) is quite general, i.e. 
it does not depend on the particular form of the functional dependence of A(t) and B{t), 
whereas the values of the parameters Y and X = {Xi, X2, X3} (or $ and X = {Xi, X2, X^}) 
are, of course, determined by this functional dependence. Using (|16D one can write the 
evolution matrix as 

U{t, 0) = P(t)e-'('^*)* , P{t + T) = P{t) . (38) 

The matrix of evolution over an integer number of periods is 

U{kT, 0) = (Ut)'' = exp[-i(o-X)A;$] . (39) 

From eqs. ([TSl) and (p6D one finds 

i3T = -i(o-X)$. (40) 

Since X is a unit vector, the eigenvalues of the matrix crX are ±1, hence the characteristic 
exponents ai^2 are purely imaginary, and the characteristic numbers ai^2 are of modulus 
one: 

ai,2 = , ai,2 = e±** . (41) 

'^The reality of these parameters is a consequence of the relation (72U{t,to)*a2 — [/(t,to) which, in turn, 
follows from the properties of the Hamiltonian H{t)* = Ti.{t) and {7i(t), 0-2} = 0. 
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(I 2) 

Let us find the normal solutions. The eigenvectors 0o of the monodromy matrix 
coincide with the eigenvectors of the matrix crX; they can be written as 





1-X.e- 



(42) 



where 



6 = arg{Xi + iXa) . (43) 

Notice that 00^"* = 0o^^ = —ia2{<f)l)^)*- The normal solutions are now found as ^(^'^^(t) = 
U{t,0)(j)i^'^\ which gives 

0(i)(t) = P(t)e-^*(*/^)0i^) , <P^^\t) = P(t)e^*(*/^)0f,') . (44) 

They form a fundamental set, which means that an arbitrary solution ip{t) can be represented 
as a linear combination of 0^^(t) and (f)^^{t) with constant coefficients: 

^(t) = P(t)[Ci e-^*(*/^)rf^ + C2 e^*(*/^Vf ^] • (45) 

The normal solutions (^) are the products of the periodic functions with periods T and 
T = (27r/$)T. Thus, the general solution (|45| ) describes modulated oscillations between the 
components of ipit) - the parametric oscillations |20|. 



Let us now discuss the parametric resonance in the system under consideration. In the 
general case of eq. (P, the parametric resonance typically corresponds to the situations when 
a solution becomes unbounded, i.e. there are some values of t for which the modulus of a 
component of ip can exceed any pre-assigned number, however large. A characteristic feature 
of the parametric resonance is that this can happen even for arbitrarily small amplitude of 
variations of the coefficients of eq. (|I|). 



In the case of the systems described by eq. (|33|), the parametric resonance corresponds 
to a situation when there exist values of t for which the modulus of a component of a solution 
can reach maximal allowed by unitarity value, unattainable in the case of the corresponding 
equation with constant coefficients. This can happen even for arbitrarily small amplitude 
of the variations of the coefficients in eq. (B3). Notice that the parametric resonance is 



undefined in the cases when eq. (|33|) with constant coefficients itself leads to maximal 
amplitude oscillations. 

To be more specific, consider the case n = 2. If the Hamiltonian of the system was con- 
stant, the Schrodinger equation ( p3D would describe the oscillations between the components 
of ijj with the frequency ujq = y/A^ + B"^ and amplitude sin^ 26*0 = B^/{A^ + B^). U 0, 
this amplitude is always less than unity. We shall now consider the case of time dependent 
A and B, but will be assuming that A{t) never vanishes during the period of evolution of 
interest 0. The transition probability in the "quasi time-independent" (adiabatic) regime 

^In the case of neutrino oscillations in vacuum, the condition A = corresponds to maximum mixing, 
while for neutrino oscillation in a matter of varying density A(t) = corresponds to the MSW resonance. 
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then never exceeds maa;{i?(t)^/(A(t)^ + 5(t)^)} < 1. The parametric resonance occurs when 
there are values of t for which the modulus of a component of ip, which was initially equal 
to zero, can reach the maximal allowed by unitarity value, corresponding to the transition 
probability equal to one. This can happen even if maxlB {tY / {A{tY + B{tY)} <ti 1 and for 
arbitrarily small amplitudes of time variations of A{t) and B{t). 

Let us find the parametric resonance condition. Assume that the initial state at t = 
is ipo = (1, 0)^. Then at t = kT we find 



cos A;$ — 1X3 sin A;$ 
ij{kT) = t/(A;T,0)Vo = I , ^ 1 • (46) 

- - 

^3 



-ie-'\ 1-X? smk<i> 



We shall now show that the parametric resonance condition is |n| 



X3 = 0. (47) 

Indeed, the transition probability reaches its maximum possible value, equal to one, when 
the survival probability |'?/'i(t)p vanishes. From (|^), in the case X3 = 0, one has ipi{t = 
kT) = cos The component ipi{t) at t = {k + 1)T is cos(/c + 1)$. It is easy to see that for 
an arbitrary nonzero value of $ there is a value of k for which cosfc$ < 0, cos(/c + 1)$ > 
or vice versa. Since all the solutions of eq. (|33D with regular coefficients are continuous, this 
means that there is a value ti, kT < ti < {k + 1)T, for which ipilti) = and the survival 
probability vanishes, i.e. the component 1^2 saturates the unitarity limit. Thus { ^7\j is the 
parametric resonance condition. 



4 Neutrino oscillations in matter with "castle wall" 
density profile 

We shall now consider applications of the Floquet theory reviewed in the preceding sections 
to neutrino oscillations in a matter of periodically modulated density. In particular, we shall 
be interested in the parametric resonance in neutrino oscillations. 

Consider oscillations in a two-flavour neutrino system. The evolution of the system in 
the flavour eigenstate basis is described by the Schrodinger equation with the Hamiltonian 
([3^ ) in which the parameters A and B are given by 

A(t) = cos 200 -^N(t), 5 = ^sin2^o, (48) 

Here G^? is the Fermi constant, E is neutrino energy, Am^ = m\ — m\, where mi 2 are the 
neutrino mass eigenvalues, and 9q is the mixing angle in vacuum. The effective density N{t) 
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depends on the type of the neutrinos taking part in the oscillations: 

{Ne for Z/e ^ V^^r 

for u^^u, 

N,-NJ2 for u.^u, ^^^^ 

-Nn/2 for l^^^r^l^s- 

Here and Nn are the electron and neutron number densities, respectively. For transitions 
between antineutrinos one should substitute —N for N in eq. (|48|) . If overall matter density 
and/or chemical composition varies along the neutrino path, the effective density depends 
on the neutrino coordinate t. The instantaneous oscillation length lm{t) and mixing angle 
6m{t) in matter are given by 

l^{t) = -K/uit) , sin 2em{t) = B/uit) , uj{t) = ^B^ + Aitf . (50) 

The MSW resonance corresponds to A{tres) = 0, sm26m{tres) = 1- 

For the parametric resonance to occur, the exact shape of the matter density profile is 
not very important; what is important is that the change in the density be synchronized in 
a certain way with the change of the oscillation phase. In particular, in p, ^ the case of the 
sinusoidal density profile was considered in which the neutrino evolution equation reduces 
to a modified Mathieu equation. In the parametric resonance was also considered for 
neutrino oscillations in a matter with a periodic step function density profile, which allows 
a very simple exact analytic solution. This solution was studied in detail in [0, Here 
we shall review this solution and its main features. 

Consider the case when the effective density N{t) (and therefore A{t)) is a periodic step 
function: 

' A^i , < t < Ti 
^N2, T,<t<T^+T2 

N{t + T)= N{t) , T = Ti + T2. (51) 

Here A^i and N2 are constants. We shall call this the "castle wall" density profile (see fig. 
2). The function A{t) is expressed by similar formula with constants Ai and A2. Thus, the 
Hamiltonian Ti.{t) is also a periodic function of time with the period T. Let us denote 

. Am2 _ G 



Nit) 



' V2 

In this notation 



V, = ^N, (^ = 1, 2). (52) 



Ai = cos2eo 6 -Vi, 5 = sin 2^0 ^, Ui = ^ (cos 2606 - Vi)^ + (sin 2606)^ . (53) 

Any instant of time in the evolution of the neutrino system belongs to one of the two 
kinds of the time intervals: 

(1) : + nT<t<T^ + nT 

(2) : Ti+nT <t <Ti+T2 + nT , n = 0, 1, 2, , ^ ^ 
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In either of these time intervals the Hamiltonian 7i is a constant matrix which we denote Tii 
and Ti.2, respectively. Let us define the evolution matrices for the intervals of time (0, Ti) 
and (Ti,Ti + T2): 

Ui = exp{-iHiTi) , U2 = exp{-iH2T2) (55) 
The monodromy matrix is then 

Ut = U2U1 . (56) 

Let us introduce the unit vectors 

ni = — (5, 0, - Ai) = (sin20i, 0, -cos2^i), 

112 = —{B, 0, - A2) = (sin 2^2, 0, - cos 2^2) , (57) 

i02 

where 6*1,2 are the mixing angles in matter at densities A^i and N2: 61 = 9m{Ni), 62 = Om{N2). 
Then one can write 

Hi = uJi{fTiii) . (58) 

Using eqs. (|55D - (^) one can obtain the monodromy matrix in the form (^) with the 
parameters Y and X given by 

Y = C1C2 - (nin2)siS2 , (59) 

X = S1C2 ni + S2C1 112 - S1S2 (ni X n2) , (60) 
where we have used the notation 

Si = sin 0j , Ci = cos ipi, (t)i = UiTi {i = 1, 2) . (61) 

Notice that the difference of the neutrino eigenenergies in a matter of density Ni is 2uji, 
so that 201 and 202 are the oscillations phases acquired over the intervals Ti and T2. The 
evolution matrix for k periods (fc=0, 1, 2,...) is given by eq. (pQ]). 

The vector X can be written in components as 

X = {(siC2 sin 26*1 + S2C1 sin 26*2) , — 3152 sin(2^^i — 26*2) , — (siC2 cos 26*1 + S2C1 cos 26*2) } . 

(62) 

Eqs. (|36|) - (|39D and (|^) - (|6^) give the exact solution of the evolution equation for any 
instant of time that is an integer multiple of the period T. In order to obtain the solution 
for kT < t < {k + 1)T one has to evolve the solution at t = kT by applying the evolution 
matrix 

Ui{t, kT) = exp[-iHi ■ {t - kT)] (63) 

for kT <t<kT + TiOT 

f/2(t, kT + Ti)Ui = exp[-^7^2 ■ {t - kT - Ti)] e^-iHiTi] (64) 



for kT + Ti <t < {k + 1)T, with Tii 2 given by eq. (|58|). 
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4.1 Parametric resonance 



Assume that the initial neutrino state at t = is a flavor eigenstate z/^. The probability 
of finding another flavor eigenstate z/^ at a time t > (transition probability) is then 
P{i^a —>■ i^bi t) = \U2i{t, 0)p. As pointed out in sec. 3, the evolution of neutrino system in a 
matter of periodically varying density has a character of parametric oscillations - modulated 
oscillations characterized by two periods, T and r = (27r/$)T. As can bee seen from eq. 



the transition probability after passing k periods of density modulation is |Tj] 



Piua ^Ub, t = kT) = (1 - Xl) sin^ $p , $p = fc$ , (65) 

where $ was defined in (0). Notice that this expression is valid for any periodic mat- 
ter density profile, irrespective of its shape. The values of X3 and of course, depend 
on this shape. For neutrino oscillations in matter with the "castle wall" density profile, 
eq. (|65|) corresponds to the evolution over an even number of layers of constant density. 
The transition probability after passing an odd number of alternating layers, which can be 
considered as k periods plus one additional layer of density Ni (the corresponding evolution 
time t = kT + Ti), is also given by eq. (p5|), the only difference being that the phase is now 



20| 



^p = k^ + ip (66) 

with 



simp = si sin29i/\l 1 — Xf , 

cos(^ = (si sin 2eiY + $2 sin 262) /^Jx."^ - X| (67) 

Eqs. (^) - ( |67D give the transition probability at the borders of the layers. 

The parametric resonance occurs when the depth of the parametric oscillations described 
by (|65D becomes equal to unity, i.e. when X3 = 0. This coincides with the condition (|47| ) 
found in sec. 3. For the case of the "castle wall" density profile under consideration it reads 
H (see eq. (H)): 

X3 = -(S1C2 cos 261 + S2C1 COS 262) = . (68) 

As follows from (^), the maximum transition probability P = 1 can be achieved at the 
borders of the layers provided that 

$p = ^ + n7r, n = 0,1,2,... (69) 

The parametric resonance condition ( |68D can be realized in two different ways. One 
possibility is that both terms on the right hand side of eq. (|68| ) vanish. This requires 
Ci = C2 = [8 -14], or 

0i = | + fc'7r, (/,2 = I + fc'V , fc',A;" = 0,1,2,... (70) 



^In refs. ^, |l^ these conditions were derived for the particular case k' = fc", which includes the most 
important principal resonance with k' = k" = 0. 
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Figure 1 : Coordinate dependence of the neutrino flavor transition probability P in a matter with 
the "castle wall" density profile, sin^ 26*0 = 0.01, 6 = lO'^^ ^y, Vi = IQ-^^ eV, V2 = 6.33 x IQ-^^ 
eV, Ti = 5.4 X 10-2, T2 = 0.1296, all distances are in units of i? = 3.23 x 10^^ eV^^ 
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Figure 2: Coordinate dependence of the matter-induced neutrino potential [Gf/V^x (density 
profile)] for the case shown in fig. 1. 
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Notice that another option, si = S2 = 0, leads to a trivial case X = 0, y = ±1, in 
which the monodromy matrix coincides (up to the sign) with the unit matrix, and the 
transition probability on the borders of the layers vanishes f\. Another possible realization 
of the parametric resonance condition is when neither of the terms on the right hand side 
of eq. (|68|) vanishes but they exactly cancel each other. 

We shall consider now the realization ( [70|) of the parametric resonance condition (|68| ) 



(the second realization will be illustrated by a numerical example in fig. 7). At the resonance, 
the transition probability for the evolution over k periods of density modulation takes a 
simple form 

P(z/„ ^u,, t = kT) = sin2[A;(2^2 - 2^1)] . (71) 

Let us first assume that the densities A^^i, A^2 are either both below the MSW resonance 
density Nmswi which is determined from GfNmsw/ = cos26'o 5, or they are both above 
it. This means that the mixing angles in matter ^12 satisfy 612 < 7r/4 or 61^2 > 7r/4, 
respectively. It is easy to see that in this case the difference 262 — 26i is always farther 
away from 7r/2 than either 26i or 26'2. Therefore in this case the transition probability 
for evolution over one period cannot exceed the maximal transition probabilities in matter 
of constant density equal to either A^i or A^2, namely, sin^ 26'i or sin^26'2. However, the 
parametric resonance does lead to an important gain. In a medium of constant density Ni 
the transition probability can never exceed sin^2^j, no matter how long the distance that 
neutrinos travel. On the contrary, in the matter with "castle wall" density profile, if the 
parametric resonance conditions (|70|) are satisfied, the transition probability can become 
large provided neutrinos travel a large enough distance. It can be seen from ([7T| ) that the 
transition probability can become quite sizeable even for small sin^ 2^i and sin^ 26'2 (i.e., in 
terms of the parameters of the Hamiltonian (0), for small max{-B^/(A(t)^ + B"^)}). This 
is illustrated in figs. 2 and 3 for the case Ni,N2 < Nmsw (the transition probability in the 
case Ni,N2 > Nmsw has a similar behavior). The number of periods neutrinos have to pass 
in order to experience a complete (or almost complete) conversion is 

It is instructive to consider the limit of small density variations, \Ni — N2\ <^ A^^i. 
In terms of the analogy with a pendulum with vertically oscillating point of support, it 
corresponds to the limit of the small amplitude of these vertical oscillations. In this limit 
9i ~ 02, and eq. (|68| ) reduces to sin(0i + ^2) = 0, or 4>i + 4>2 = kn. Since (pi = UiTi, this 
condition can be written as 

n = — = —, where u = ^i— + 1^2^ ■ (73) 



^We do not consider the trivial cases of the MSW resonance for which X3 = because cos 26i = and 
Si = ±1, z = 1 or 2, or cos 20i ~ cos2^?2 ~ 0. These cases correspond to A{tres) ~ 0. 
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Figure 3: Same as in fig. 1 but for 6 = 10"^^ eV, V2 = 10"^^ eV, Ti = 5.4 x 10"^ T2 = 5.4 x 10 
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Figure 4: Coordinate dependence of the matter-induced neutrino potential for the case shown 
in fig. 3. 
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Figure 5: Same as in fig. 3 but for V2 = Vi {V{t) = V\ = const). 



This coincides with the famihar parametric resonance condition in the case of the small- 
amphtude variations of the parameter of the system, eq. (|3^). It is important to notice that 



the condition (|73| ) does not depend on the amphtude of the density modulation, A^^i — N2. 
This illustrates the point that we emphasized in sec. 3 - the parametric resonance can occur 
even for arbitrarily small amplitude of the variations of the parameters of the system. Of 
course, the smaller this amplitude, the longer the evolution time for the total conversion. 

Consider now the case Ni < Nmsw < ^2 (^1 < ti / 4: < 62) ■ The transition probability 
over n periods at the parametric resonance is again given by eq. (0). However in this case, 
for 02 > tt/ 4 + 61/2 (which is always satisfied for small mixing in matter), one has sin^(2^2 — 
29i) > sin^20i, sin^2^2- This means that even for the time interval equal to one period 
of matter density modulation the transition probability exceeds the maximal probabilities of 
oscillations in matter of constant densities Ni and N2- The case A^i < Nmsw < ^2 is 
illustrated in figs. 3, 6 and 7. 

Figs. 3 and 4 show the importance of the phase relationships in the case of the parametric 
resonance. In these figures the coordinate dependence of the transition probability and 
matter density profile are shown for a specific case in which conditions ( [701 ) are fulfilled. It 
can be seen from these figures that the probability increase during the time intervals T2, 
which correspond to the effective matter density N2, is very small, and, in addition, in this 
case T2 -C Ti. One could therefore conclude that the evolution during these very narrow 
intervals is unimportant. However, this conclusion is wrong: if one removes the "spikes" in 
the matter density profile of fig. 4, i.e. replaces it by the profile N{t) = Ni = const, the 
resulting transition probability will be very small at all times (fig. 5). 

Further examples of the parametrically enhanced neutrino oscillations can be found in 
figs. 6 and 7. Figures 1, 3 and 6 correspond to the reahzation (^) (ci = C2 = 0) of the 
parametric resonance condition (^); fig. 7 illustrates the realization in which the two terms 
in X3 cancel each other. 
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Figure 6: Solid curve: Coordinate dependence of the transition probability P for the case of 
total conversion over 5 periods of density modulation (10 layers). Dashed curve: the same for 
the case of total conversion over 3 layers. The kinks correspond to the borders of the layers of 
different densities. The curves were plotted for the realization (|70| ) (ci = C2 = 0) of the parametric 
resonance condition. 
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5 Parametric resonance in neutrino oscillations in the 
earth 



5.1 Evolution of oscillating neutrinos in the earth 

The earth consists of two main structures, the mantle and the core, which can to a very 
good approximation be considered as layers of constant density. We shall consider neutrino 
oscillations in the earth in this two-layer approximation. Neutrinos coming to the detector 
from the lower hemisphere of the earth at zenith angles G in the range cos© = (— 1) -^ 
(—0.837) (nadir angle G„ = 180° — < 33.17°) traverse the earth's mantle, core and then 
again mantle, i.e. three layers of constant density with the third layer being identical to 
the first one. Therefore such neutrinos experience a periodic "castle wall" potential, and 
their oscillations can be parametrically enhanced. Although the neutrinos propagate in 
this case only through three layers ("1.5 periods" of density modulation), the parametric 
enhancement of the transition probability can be very strong. 

The evolution matrix in this case is U = U1U2U1. It can be parametrized in a form 
similar to that in eq. (pGf): 

U = Z - ia-W , Z^ + W^ = l. (74) 

The matrix U describes the evolution of an arbitrary initial state and therefore contains 
all the information relevant for neutrino oscillations. In particular, the probabilities of the 
neutrino flavor oscillations P and of z/2 ^ oscillations (relevant for the oscillations of 
solar and supernova neutrinos inside the earth) are given by 



P = W^ + Wi , Pae = sin^ ^0 + Wi{Wi cos 29^ + sin 2^o) • (75) 
Equivalently, the probability P can be described by eqs. (p5| ) - ( |67D with k = 1. 



We have now to identify the effective densities A^^i and N2 with the average matter 
densities Nm and Nc in the earth's mantle and core, respectively; similarly, we change the 
notation Vi^2 Vrnfi, ^1,2 ^ 0m,c and 6*1^2 ^ ^m,c- 

In the two-layer approximation, the parameters Z, W have a very simple form [IIJ 



Z = 2 COS0m ^ — COS0C , (76) 

W = (2 sin 0m sin 2Qra Y + sin 0c sin 26c , , — (2 sin 0^ cos 26m Y + sin 0c cos 2^c)) • (77) 
Here the vector W was written in components; the parameter Y was defined in ([59|). If the 



parametric resonance condition (|68D is satisfied through the realization ([TOD, the neutrino 
flavor transition probability takes the value [111], O 



P = sin^(2^c-4^m), (75 



^Different but equivalent expressions can also be found in p3, [IJ 
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whereas the probabihty of the z/2 ^ Ve transitions is 



m 



P2e = sm\29, - Adm + ^0) • (79) 

These probabihties can be close to unity (the arguments of the sines close to vr/2) even 
if the amplitudes of neutrino oscillations in the mantle, sin^ 2^^, and in the core, sin^ 26c, 
are rather small. This can happen if the neutrino energy lies in the range < E < E^, 
where Em and Ec are the values of the energy that correspond to the MSW resonance in the 
mantle and in the core of the earth. This condition is equivalent to A^^ < Nmsw < ^c- In 
the case of small mixing angle MSW solution of the solar neutrino problem, sin^ 29o < 10~^, 
and practically coincides with P unless both probabilities are very small. 

The trajectories of neutrinos traversing the earth are determined by their nadir angle 
Qn = 180° — G. The distances Rm and Rc that neutrinos travel in the mantle (each layer) 
and in the core are given by 

Rm = R (cos e„ - ^rVi?2 - sin^ e„ ^ , R^ = 2R y'rVi?2 - sin^ e„ . (80) 

Here R = 6371 km is the earth's radius and r = 3486 km is the radius of the core. The 
matter density in the mantle of the earth ranges from 2.7 g/cm^ at the surface to 5.5 g/cm? 
at the bottom, and that in the core ranges from 9.9 to 12.5 gjcvn? (see, e.g., [^). The 
electron number fraction Y^. is close to 1/2 both in the mantle and in the core. Taking the 
average matter densities in the mantle and core to be 4.5 and 11.5 gjav? respectively, one 
finds for the v^. ^ i/^ oscillations involving only active neutrinos the following values of 
Vra and Vc- Vm = 8.58 x 10~^^ eV, Vc = 2.19 x 10~^^ eV. For transitions involving sterile 
neutrinos Ue ^ Vg and v^^^r ^ t^s, these parameters are a factor of two smaller. 



5.2 Parametric resonance conditions for neutrino oscillations in 
the earth 

If the parametric resonance conditions (|70D are satisfied, strong parametric enhancement 
of the oscillations of core crossing neutrinos in the earth can occur |]11], O, IT3|, IT^, |T^, ff6 



see fig. 8. In some cases, condition ( p9D can also be fulfilled, and the parametric resonance 
leads to a complete flavour conversion for neutrinos traversing the earth. 



We shall now discuss the resonance conditions ([701). The phases 0^ and 0c depend on 
the neutrino parameters Am^, 6q and E and also on the distances Rm and Rc that the 
neutrinos travel in the mantle and in the core. The path lengths Rm and Rc vary with the 
nadir angle; however, as can be seen from (pOD , their changes are correlated and they cannot 
take arbitrary values. Therefore if for some values of the neutrino parameters a value of the 
nadir angle 6„ exists for which, for example, the first condition in eq. (^) is satisfied, it is 
not obvious if at the same value of 6„ the second condition will be satisfied as well. In other 
words, it is not clear if the realization (^) of the parametric resonance condition ( |68D is 
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Figure 8: Solid curve: transition probability P for <^ v^^r oscillations in the earth as a function 
of the distance t (measured in units of the earth's radius R) along the neutrino trajectory. 5 = 
/^ni^/AE = 1.8 X lO'^^ eV, sin^ 2(9o = 0.01, e„ = 11.5°. Dashed curve: the same for a hypothetical 
case of neutrino propagation over full two periods of density modulation {tmax = "^{Rm + Rc))- 



possible for neutrino oscillations in the earth for at least one set of the neutrino parameters 
Am^, ^0 and E. However, as was shown in [^, not only the parametric resonance 
conditions ( [70D are satisfied (or approximately satisfied) for a rather wide range of the nadir 
angles covering the earth's core, they are fulfilled for the ranges of neutrino parameters 
which are of interest for the neutrino oscillations solutions of the solar and atmospheric 
neutrino problems. In particular, the conditions for the principal resonance {k' = k" = 0) 
are satisfied to a good accuracy for sin^ 29o ^ 0.1, 6 ~ (1.1 -^ 1.9) x 10~^^ eV^, which includes 
the ranges relevant for the small mixing angle MSW solution of the solar neutrino problem 
and for the subdominant z/^ u^, and ^ Vt oscillations of atmospheric neutrinos. 

The fact that the parametric resonance conditions (^) can be satisfied so well for 
neutrino oscillations in the earth is rather surprising. It is a consequence of a number of 
remarkable numerical coincidences. It has been known for some time that the potentials Vm 
and Vc corresponding to the matter densities in the mantle and core, the inverse radius of the 
earth , and typical values of (5 = /\m? /AE of interest for solar and atmospheric neutrinos, 
are all of the same order of magnitude - (3 x 10~^^ - 3 x 10"^^) eV (see, e.g., [^, |lT], ^). It 
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is this surprising coincidence that makes appreciable earth matter effects on the oscillations 
of solar and atmospheric neutrinos possible. However, for the parametric resonance to take 
place, a coincidence by an order of magnitude is not sufficient: the conditions (|70D have to be 



satisfied at least within a 50% accuracy |T^. This is exactly what takes place. In addition, 
in a wide range of the nadir angles B^, with changing 0^ the value of the parameter 6 
at which the resonance conditions (|70D are satisfied slightly changes, but the fulfillment of 
these conditions is not destroyed. 

Even more surprising is the fact that the second realization of the parametric resonance 
condition (|68l), the one in which the two terms in X3 cancel each other, is also possible for 



neutrino oscillations in the earth . This requires a very subtle tuning between the values 
of neutrino parameters ^0 and S = Am'^/AE on one hand and the effective matter densities 
Nm and Nc and neutrino pathlengths Rm and Rc in the mantle and in the core on the other. 
This looks very contrived, and yet turns out to be possible. Moreover, for a number of the 
values of the neutrino parameters, the condition ( p9| ) can also be fulfilled and a complete 
flavour conversion for neutrinos crossing the earth is is possible [|19|. 



6 Discussion and conclusion 

We have reviewed the Floquet theory of linear differential equations with periodic coef- 
ficients and discussed its applications for neutrino oscillations in matter of periodically 
modulated density. In particular, in the case of two-flavour oscillations, we have shown that 
the evolution of the system takes a form of parametric oscillations - modulated oscillations 
characterized by two periods, the period T of density modulations and r = (27r/$)T, where 
±i$ are the characteristic exponents. We have also discussed the parametric resonance in 
neutrino oscillations and have shown that, irrespective of the shape of the periodic density 
modulations, the parametric resonance condition is X3 = 0, where is one of the pa- 
rameters that determine the monodromy matrix given in eq. (|36|) . We have also reviewed 
the exact solution in the case of 2-flavour neutrino oscillations in a matter of periodic step- 
function density profile, obtained in P, [1^, and discussed its relations with the Floquet 
theory. This solutions allows one to find explicit formulas for the parameters entering into 
the monodromy matrix, including the characteristic exponents and X3. We discussed the 
implications of this exact solution for the oscillations of neutrinos inside the earth, the den- 
sity profile of which can to a very good accuracy be approximated by a piece of the periodic 
step-function profile. We concentrated on possible parametric resonance effects in neutrino 
oscillations in the earth. 

Besides being an interesting physical phenomenon, the parametric resonance in neutrino 
oscillations can provide us with an important additional information about neutrino prop- 
erties. Therefore experimental observation of this effect would be of considerable interest. 
The prospects for experimental observation of the parametric resonance in oscillations of 
solar and atmospheric neutrinos traversing the earth were discussed in [0, To a large 
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extent, these prospects depend on the values of some of neutrino parameters which are not 
well known yet. The bottom line is that such an experimental observation is difficult but 
may be possible. 

Parametric enhancement may lead to noticeable effects in oscillations of supernova neu- 
trinos in the earth, resulting in characteristic distortions of the spectra of neutrinos crossing 



the earth's core ||2^. A sufficiently accurate measurement of the supernova neutrino spec- 
trum would, however, require a relatively close supernova (L < 10 kpc). 

An interesting possibility to study the parametric effects in neutrino oscillations would be 
very long baseline experiments with intense neutrino beams produced at neutrino factories 
(for discussions of neutrino oscillations experiments at neutrino factories see, e.g., |^). For 
baselines larger than approximately 10700 km, neutrinos would cross the earth's core, and it 
would be possible to probe the parametric resonance effects in v^. v^(t) oscillations. Notice, 
however, that the present feasibility studies concentrate on relatively short baselines, a few 
thousand km In addition, the expected average energies of neutrino beams [E > 20 



GeV) in the currently discussed experiments are somewhat higher than what would be 
desirable in order to study the parametric effects. 

As we have seen, observing the parametric resonance in oscillations of solar, atmo- 
spheric or supernova neutrinos in the earth or in experiments with neutrino factories is not 
an easy task. Can one create the necessary matter density profile and observe the paramet- 
ric resonance in neutrino oscillations in the laboratory (i.e. short-baseline) experiments? 
Unfortunately, the answer to this question seems to be negative: this would require either 
too long a baseline or neutrino propagation in a matter of too high a density (see [T^, IE 



for details). One can conclude that the sole presently known object where the parametric 
resonance in neutrino oscillations can take place is our planet, as was first pointed out in 



11, 12 
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